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Abstract 

Interrelation between Thorn's catastrophes and differential equations revisited. 
It is shown that versal deformations of critical points for singularities of A,D,E type 
are described by the systems of Hamilton-Jacobi type equations. For particular 
nonversal unfoldings the corresponding equations are equivalent to the integrable 
two-component hydrodynamic type systems like classical shallow water equation, 
dispersionless Toda system and others. Peculiarity of such integrable systems is 
that the generating functions for the corresponding hierarchies, which obey Euler- 
Poisson-Darboux equation, contain information about normal forms of higher order 
and higher corank singularities. 
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1 Introduction 



Connection between catastrophe theory ( see e.g. [1, 2, 3, 4, 5]) and theory of singularities 
for differential equations has been studied for more than fourty years. It was R.Thom 
himself [1] to observe that the singularities developed by solutions of the wave equation 
are related to the unfolding of singularities of functions. In 1968 E. Zeeman noted [6] 
that the breaking waves phenomenon (gradient catastrophe) is associated with the hyper- 
bolic umbilic catastrophe. Within the formalism of Lagrangian submanifolds the above 
interrelation has been studied by Guckenheimer [7], Arnold [8, 9] and other authors ( 
see e.g. [10, 11, 12, 1-3, 14, 15, 16] and references therein). Some aspects of the tangent 
fields for singularity unfoldings in connection with deformations of complex spaces have 
been considered in [17] . In topological quantum field theory [18] unfoldings of singu- 
larities and corresponding Witten-Dijkgraaf-Verlinde-Verlinde (WDVV) equations have 
appeared within the study of the deformed Ginzburg-Landau models [19, 20]. A connec- 
tion of the topological field theory and Saito's approach in singularity theory [21] with 
Frobenius manifolds [22]and associated WDVV and Hirora equations has been discussed 
in [22, 23, 24, 25, 26, 27, 28]. It is well established nowadays that the base spaces of versal 
unfoldings of singularities may carry rich family of algebraic and geometrical structures 
including various types of differential equations. 

Within different studies an appearance of normal forms of singularities has been ob- 
served recently in the papers [29, 30, 31, 32, 33, 34]. In [30, 31, 32] cusp and umbilic 
catastrophes have arozen within the study of the behaviour of the full dispersive integrable 
equations near the points of gradient catastrophe. In contrast, an analysis performed in 
[29, 33, 34] deals directly with the dispersionless integrable systems of certain class. In 
the papers [33, 34, 35] it was shown that the hodograph equations for these systems are 
nothing else than the equations defining critical points of the functions W which obey 
Euler-Poisson-Darboux equations. Such functions W have several important properties: 
they are the generating functions for the whole integrable hierarchies, they drastically 
simplify an analysis of the singular sectors for these hierarchies. But, perhaps, more in- 
teresting fact is that they contain information about the normal forms of singularities of 
functions from the catastrophe theory. The observations made in the papers [29, 33, 34] 
deal with the singularities of An type and umbilic type. 

In the present paper we will discuss one particular aspect of the relation between 
the singularity theory and differential equations. Let the function F{x;t) = Fo{x) + 
'^2=o^k{x)tk with certain functions efc(x) and deformation parameters tk defines the un- 
folding of the singularity with the normal form given by the function Fo{x) of m variables 
Xi, Xm.. We will be interested in differential equations which govern the dependence of 
the critical points Ui,i = 1, ...,m for the function F{x]t) on deformation parameters. We 
will show, using an elementary technique and addressing mainly to nonexperts, that 

1. For the versal deformations of the A,D,E singularities ( m = 1 or 2) the dependence 
of their critical points on outside the catastrophe sets is described by the systems of 
mn differential equations of the first order. These systems imply that Ui = ip^. and the 
function ip obeys the systems of Hamilton- Jacobi type equations. For An singularities 
these equations are of the form 
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V'i, = (<^ij, A; = 2,...,n-1 

or 

= {u'')t^, k = 2, ...,n - 1, 

which are the Burgers-Hopf and n-3 higher Burgers-Hopf equations. For D„, Eq, E^, Eg 
cases (m=2) the corresponding equations are 

V'i, = ez(Vfi,V't2), / = 3,...,/i-l 

or 

«, = ^^, -U = ^^^ / = 3,...,;.-l 

where u = ui,v = U2 and ei{x,y) are elements of a basis of the local algebra Qfq for 
critical points of the dimension In all these cases for the versal deformations 

F{u; dt) = dip. (1) 

2. For the particular nonversal unfolding of the umbilic singularities, for which 
F{x, y; t) = Fo{x, y) + t^F^x +tix and the constraint {F^^ — 6Fyy){u, v) = with 5 = ±1 
is verified, the critical points u{t) and v{t) obey integrable systems of hydro dynamic type 
of the form 

'"^ DC) (2) 



V I ^ \ Sv, au J \ V , ^ 

/ t:j \ / \ / ti 

For certain values of a this system represents the well-known integrable systems like the 
classical shallow waver equation [a = 2) and dispersionless Toda system (a = 0). In this 
case the catastrophe set defined by the equation A -^ {F^xFyy — Fxy){u, v) = decomposes 
in two components. For such integrable systems the corresponding hodograph equations 
coincide with the equations for critical points u and v of the functions W obeying the 
Euler-Poisson-Darboux equations. For this unfolding and other nonversal unfoldings of 
such type the relation (1) is valid too. 

3. It is observed that the generating functions W{x;t) = J2k>i^kWk{x) for the inte- 
grable n-component hydrodynamic type systems contain information about normal forms 
of singularities of higher orders and coranks. It is shown that in the two-component case, 
namely, for the system (2) 

W3^x^± 3xy^, W4~x^ + dx^y^ + y\ 

i.e. the functions W3 and W4 are proportional to the normal forms of and Xg singu- 
larities of corank two, respectively. In the three-component case 

x^ + + + dxyz, 
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that is the normal form of the unimodular singularity Pg (or T3 3 3) of corank three. De- 
pendence of the critical points for this singularity on deformation parameters is governed 
by the dispersionless three-component coupled KdV equation. 



2 Thorn's catastrophes 

First, we recall briefly the basic facts about unfolding of singularities (see [1, 2, 3, 4, 
5]). Let the normal form of the singularity of the corank m is given by the function 
Fq{xi, ...,Xm) ■ An unfolding ( deformation) F{x;X) of the function Fq is a germ of a 
smooth function F : i?"' x i?" -> at the point (0,0) such that F{x;0) = Fo{x). The 
space i?" of the second argument of F is called the base space of unfolding and Ai, A„ are 
called the parameters of unfolding. Dimension of the base space for different deformations 
of the same function can be different. Critical set of F{x; A) is given by the solutions of 
the equations Fx-{u; \) ^(m;A) = 0,i = l,...,m and the catastophe set is defined 
by vanishing of the Hessian (w) = 0. The tangent space of the orbit of the germ 

Fq{x) is its gradient ( Jacobian) ideal Ij = ( -Fqxi, -Foa;,„) which consists of all germs of 
the functions of the form "^^jZ^ hi{x)FQx.. Local algebra Qfo of the critical point is the 
quotient R[x]/Ij. It's dimension fi is called the Milnor number ( or multiplicity) of the 
critical point. If the functions cq = l,ei{x), e^_i(x) form a basis in Qfo^ then a versal 
( R- versal) unfolding of Fq{x) can be represented in the form 

F{x-t) = Fo{x) + Y,ek{x)tk. (3) 

fc=0 

Here are the deformation parameters and hence the initial velocities of deformation 
Vk -T- Ft^{x; t)\t=o = ek{x), k = 0,1, /i — 1. Nonversal unfoldings of Fq{x) may depend 
on less than /i parameters or have other forms. 

Let us begin with the celebrated seven Thom's catastrophes. The corresponding 
functions Fq and their versal unfoldings are (see [1, 2, 3, 4, 5]): for corank one An type 
catastrophes ( fold, cusp, swallow tail, butterfly) 

n-l 

F„^(x;^)=x"+l + ^^fca;^ n = 2,3,4,5, (4) 

fc=0 

for corank two hyperbolic and elliptic umbilic catastrophes ( D^) 

F^{x, y; t)=x^± 3xy^ + ts{x^ T y^) + + hx + to, (5) 
and for the parabolic umbilic catastrophe -D5 

F5(x, y; t) = x^y + + t^y"^ + ^x^ + t2y + hx + to- (6) 

We begin with the An type singularities. Critical points u of the functions (4) are 
defined by the equation 
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n-1 

Ft - Ft{u- t) = {n + l)u^ + ^ khu'-' = (7) 

k=l 

where etc denotes the derivative w.r.t. x. Calculating the differential of (7), one gets 

n— 1 n—1 

F^Ju + ku'-'dtk = Y^iF^^^Ut, + ku'-')dh = 0. (8) 

k=l k=l 

Hence, the critical points u{t) obey the system of equations 

Ftu^^t, + ku'-' = 0, /c = 1,2,..., (9) 
Outside the catastrophe sets where -F^u^ 7^ one has the system 

= -^^' fc = l,2,...,n-l. (10) 

nuu 

On the catastrophe sets defined by the condition -F^u„ = the velocities Mf^. of variations 
of the positions of the critical point becomes unbounded which is a standard manifestation 
of the catasrophe. 

For n = 2 one has a single equation while at n > 3 the system (10) is equivalent to 
the system 

«fi = -t;^, ut^ = {u'')ti, A; = 2,3, ...,n-l (11) 

nuu 

Equations (11) imply that there exists a function (p such that u = if^^ and that the n-2 
last equations (11) take the form 

^t, = M\ /c = 2,3,...,n-l. (12) 

As a consequence, evaluating the infinitesimal variation of the function (4) on the critical 
set, one obtains 

6F:^{u) ~ F^{u- t + dt)- F^{u- 1) = F^{u- dt) = d^. (13) 

Equation (11) with /c = 2 is the well-known Burgers- Hopf (BH) equation (see e.g. [36] 
) while equations (11) at A; > 3 represent n — 3 higher flows commuting with them. In the 
form (12) this system is the family of Hamilton- Jacobi type equations. The fact that the 
deformations of the critical points for the first four Thom's catastrophes are governed by 
the BH and higher BH equations has been discussed in [29] and [26]. 

Umbihc type singularities can be considered simultaneously. Indeed, let us take the 
function 

F(x, y; t) = Foix, y; t) + Uy^ + t^x^ + t^y + hx + to- (14) 
For singularity Fq = x^ + Sxy"^ and ^4 + ^3 = 0. For singularity Fq = x^ — Sxy"^ 

and ^4 — ^3 = and in case Fq = x'^y + y^. Critical set is defined by the equations 
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F^{u, v; t) = Fq^{u, v) + 2t3X + ti = 0, 



(15) 



Fy{u, V- 1) = Foy{u, v) + 2Uy + = 0. (16) 

Calculating differentials of these equations, taking into account that in the critical set 
du = Ylk=i'^tk^^k, dv = Ylk=i'^tk^^k and solving systems of linear equations, one gets 

F F 
wti = = (17) 

Fuv -^UU / 1 o \ 

Ut2 = vt2 = — (18) 

Ut3 = -2m^, = 2u— , (19) 

uu = 2v^, v,, = -2v^ (20) 

where F^u = F^xiu.v) etc and A = F^uFw — {FuvY- These equations govern the 
dependence of the critical points on deformation parameters. On the catastrophe set 
A = derivatives of u and v become infinite and one observes a very fast ( catastrophic) 
change of positions and number of critical points. 

The system of equations (17)-(20) can be represented in different equivalent forms. 
For example, it implies that 

Ut2=vti, (21) 

and 

= {u^)tu = (22) 

uu = {v\,, vt, = {v\^. (23) 

Hence, there exists a function tp such that u = ipf-_^,v = ip^^ and equations (22) and (23) 
are reduced to two equations 

^t, = i¥>tj\ ^u = i¥>tf- (24) 

It is easy to see that the system (17)-(20) is equivalent to the following system for the 
function ip 

Pvv Fuv F^u fc\r\ 

Vhh = Vht, = 'ft2t2 = — ^> (25) 

'^ = ^t„ v = ipt^, (26) 
^t, = i'fj', ^u = {^t,?- (27) 
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As the consequence of equations (26), (27) one has 

fc=4 

5F{u, v) ^ F{u, V] t + dt)- F{u, v; t) = F{u, v; dt) = ^ (ftjtk = dip (28) 

fc=0 

with = 1. 

Equations (17)- (27) describe dependence of critical points of singularity on defor- 
mation parameters. To get the corresponding equations for singularities one can pass 
in equations (17)-(27) to new variables t± = ^(^4 ± ts) and then impose the constraints 
t+ = or t_ = or use directly the relation (28). Indeed, the functions F^ for 
singularities are F"^ = F{x,y;t)\t^=zft3- So, 

k=3 

dip^ = J2 'Ptdtk = d^l^^^,^ = {{^if T {^if)dh + ^idh + ipidU + dto. (29) 

A;=0 

Hence, for the singularities one has the equation 

^t = {^i?T{^i?. (30) 

plus the corresponding equations (25), (26). In terms of m = and v = ^pf^ equation 
(30) assumes the form of the 2-|-l-dimensional hydrodynamical type system 

ut-s = (n^=Ff^)ti, 

= iu'Tv')t,. (31) 

The catastrophe sets A = for the umbilic singularities discussed above are those 
subsets of variables tk for which solutions of considered systems of differential equations 
exhibit gradient catastrophe Ut^,Vt^. 00. 



3 A,D,E singularities 

In his seminal paper [S] Arnold proved that the list of functions with simple (without 
moduli) degenerate critical points consists of two infinite series A„(n > 1), Dn{n > 4) and 
three exceptional cases EQ,Ef,Es. The corresponding normal forms and versal unfoldings 
are as follows: 

An-. F(x;t) = a;"+^ + t„_ia;"-^ + t„_2x"-2 + ... + tiX + to, (32) 
Dn-. F{x,y;t) = x^y±y''-^ + tn^iy''-^ + ...t2y + tix + to, (33) 
Ee: F{x,y;t) = x^±y^ + t5xy^ + Uxy + t3y'^ + t2y + tix + to, (34) 
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E^ : y- 1) = + xy^ + hxy + tg'/ + Uy^ + ^32/^ + + tia; + to, (35) 



F{x, y; t) = x^ + y^ + tjxy^ + t^xy^ + tgxy + Uy^ + tg?/^ + t2y + hx + to- (36) 



Equations governing the dependence of the critical points of An type has been found, 
in fact, in the previous section. It is sufficient to extend the formulas (7)-(13) to arbitrary 
n. Thus, in this case (for n > 3) one has the system of equations 



= 2, 3, n — 1 



(37) 



or 



V.,^^ = (y.,J^ A; = 2,3,...,n-1 (38) 

and F{u; dt) = d(p. The system (37) contains the BH equation together with its — 3 
higher flows. Thus, the system of differential equations governing versal deformations of 
the critical points for the entire An series with unbounded n represents the whole infinite 
BH hierarchy (see [29] and [26]). 

In order to consider D and E cases let us write the functions (33)-(36) as 



F{x, y; t) = Fo(x, y) + ^ ek{x, y)tk 



k=0 

with corresponding Fq, Ck and yU. Critical points are defined by the system 



Fu — Fqu + Gkutk — 0, 



k=0 
M-1 



Fy — Fq^ + ^ Ckvtk — 0. 

k=0 

Differentiating these equations w.r.t. ti, one gets the system 

/ = 1, yU — 1 



Fuu Fuv 
Fuv Fvv 



and consequently 



1 

' A 



= --r 



eiu 
eiv 



Fuv (^Iv 



I = 1, fj, — 1. 



From equations (42) it follows that 



dei{u,v) 



dem{u,v) 



eiuUtr, 



dtm dti 
Hence, there exists a function ip{t) such that 



eivVtr. 



(39) 



(40) 
(41) 



(42) 



(43) 



(44) 
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ei{u,v) = ipt^, / = 1, - 1. (45) 

For D and E cases one has Cq = 1, ei = x, 62 = y. Thus, ^9^^ = l,u = ipf_^,v = ip^.^ and 
equations (45) are of the form 

^u= (^lift^^ft^)^ / = 3,...,/i- 1. (46) 
In terms of u and v one has the systems 



dei{u,v) dei{u,v) 
n, = ^^, -, = ^^^ / = 3, ...,/.-!. (47) 

In virtue of (44) the systems (46), (47) describe commuting flows. From the above equa- 
tions one also concludes that for all D and E cases 

5F{u,v) ^ F{u,v;dt) = dip. (48) 

The relations of the type (48) have appeared earlier within the study of semiversal un- 
foldings of hyper surf ace singularities ( see e.g. [25], Chapter 5). 

The systems (46) or (47) governs motion of critical points for versal deformations of 
D and E singularities. 

Concretely, for the Dn case one has the system 

¥^t, = iVt,)'-\ fc = 3,4,...,n-l (49) 

and u = ip^^,v = (p^^ . Effectively, it is again the family of the BH and higher BH equations 
as in An case. 

For Eq case one has the equations 

or three systems 



For Ej singularity (35) one has 



(51) 



or equivalently 



Finally, for Eg singularity one has the system of five equations 



(53) 



fh = (<^t2)^ fu = (<^t2)^ ft, = ft^ft2^ fte = <^^l(<^^2)^ ft^ = ftAft2f (54) 
or the systems 
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= )t2, ' = {v^)t2, ' vts = {uv)t^^ ' vt^ = {uv^)t2, ' vtr = {uv^)t2, 

Interrelations between the systems (50)-(55) and their 1+1- dimensional reductions are of 
interest. It would be also of interest to analyse a connection between apparently different 
systems of equations corresponding to different choices of a basis for the local algebra Q Fq 
of the critical points. 

Analysis presented above can be easily extended to the singularities and their versal 
unfoldings of any corank m. One derives that the set of critical points Ui{i = 1, m) are 
components of the gradient Ui = ^PtM = l,..-,m) and the function ip{ti, ...,t^_i) obeys 
the system of equations 

ft,=ek{Vt,^--^'PtJ, k = m + l,...,fi-l (56) 

where efc(xi, Xm), k = 0, 1,...,// — 1 form a basis of the local algebra Qfo ioi the critical 
point of the function -Fo(^i) ^n)- Equivalently, one has the systems 

du, _ dekiu^,...,u^) k-m + 1 a-1 (57) 

and also 

F{ui, ...,Um]dt) = dip. (58) 
Properties of these equations will be discussed elsewhere. 

4 Nonversal unfoldings and integrable systems of hy- 
drodynamic type. 



Construction given in the previous section is applicable to nonversal unfoldings too. For 
nonversal unfoldings, even if the infinitesimal deformation is the form (39), the functions 
Cfc may not form a basis of the local algebra or number of parameters of deformation 
may be less than /i. Nevertheless, it is not difficult to show that for such unfoldings of 
corank two singularities the general formulas (40)-(45) remain unaltered, only number of 
the variables can be different. 

Here we will consider a particular unfolding of the umbilic singularities given by the 
function 

F{x, y- ti, ta) = ax^ + Sxy^ + t3{(3x'^ + -fy^) + hx (59) 

where a, /3, 7 are parameters. At a = ±1 and /9 ± 7 = this unfolding is equivalent to 
the versal unfoldings of singularity with the "frozen" parameters ^2 = and = 0. 
At a = the function (59) represents the two-dimensional part of the infinite dimensional 
deformation of the germ Sxy"^ of the critical point with infinite multiplicity [8]. 
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We will be interested in the subclass of unfoldings (59) for which 63 = + 7?/^ 
belongs to the tangent space of the germ ax^ + Sxy"^, i.e. 63 ~ Fq^. This contraint is 
verified if = (3. In this case the function F after trivial rescaling takes the form 

F{x, y; ti, ts) = ax^ + 3xy'^ + ts{ax'^ + y"^) + tix. (60) 
Repeating the calculation performed in the previous section, one gets the system 

F F 

^ vv ^ uv //^ 1 \ 

Ut, = — ^, Vt, = (61) 

A A ' ^g2) 

Using (61), one can rewrite the last system as 

Ut, = {au^ + v%^, 

vt, = -2v^ + 2auvt^. 

This system is not of hydrodynamical type. In absence of ^2 equations (18) we cannot 
express in terms of u and v and their derivaties. A way to get the hydrodynamical 
type system is to impose a constraint on F. It is not difficult to show that among the 
constraints having form of linear relation between second order derivatives of F only the 
constraint = 5F^^ where 5 is a constant is admissable ( nontrivial). Under this 
constraint = 6^ = —5ut, and the system (63) becomes 

M = 2 f M f M . (64) 

The first equation of this system, i.e. Ut, = {au"^ + v'^)t, implies the existence of a function 
ip{ti,t3) such that u = ip^_^ and au"^ + v"^ = ip^^. As the consequence one has 

F(n, V] dt) = F{u, v;t + dt) — F{u, v; t) = t^i^au^ + v'^) + tiu = dip. 

Among two parameters in this system only one is relevant. Indeed for nonvanishing 
S the rescalling u ^ u,v'^ ^ — > — converts the system (64) into the 

same system with 6 = ±1 leaving only the parameter a free. For nonvanishing a one can 
convert the system (64) into that with a = 1 and free 6. Since we will consider the case 
a = among the others the first choice is preferable. 

The system (64) with S = ±1 and particular values of a coincides with some well- 
known integrable hydrodynamic type systems. Indeed, at a = one has the system 



Ut, = 2vvt^, Vt, = 2Svut, (65) 

which is equivalent to the dispersionless Toda equation $43^3 = 4(5(exp ^)titi for $ = log v^. 
For a = 2 and 5 = 1 it is the hyperbolic one-layer Benney system ( classical shallow water 
equation) for variables u and v"^ ( see e.g. [36] and [37]). At a = 2 and 5 = —1 it is the 
elliptic one-layer Benney system which is equivalent to the dispersionless Da Rios system. 
For a = 6 = ±1 the corresponding system Ut, = (±m^ -|- w^)ti, Vt, = ±2{uv)ti decomposes 
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into two BH equations {u ± ^/±v)^ = =f((w ± "s/i"^)^)*!- other values of a the system 
(64) is equivalent to the system ( t = — 2at3, x = ti, p = v°'^) 

Px 

Pt + {pu)x = 0, ut + UxU + — = 0, (66) 

P 

which describes one-dimensional motions of ideal barotropic gas with the density p and 
pressure p = -^^^P (see e.g. [ib\). 

In terms of Riemann invariants (5^ = u + -^v, (32 = u — -^v the system (64) is of the 
form 

/3u3 = 1^1 + ^2) + P2t, = 1^1 {e{Pi + ^2) + /32)/32*„ (67) 

where ^ = |(|f| ~ !)• This form of the system (64) clearly shows that only the ratio is 
relevant. We will put 5 = 1 in what follows. The hodograph equations are of the form 
Fp^ = 0, Fp^ = while the constraint = F^^ becomes -F/Sj/Sj = 0. Under this condition 
A = Fp^p_^Fi3^p^ and hence the catastrophe set where the solutions of (64) or (67) exhibit 
gradient catastrophe decomposes into two pieces Ff^^p^ = and Fjs^fs^ = 0. 

The system (67) represents a two-component case of the so-called e-systems discussed 
in [39] . In a different manner the systems (64)- (66) arozen within the study of hydrody- 
namic type systems associated with the two dimensional Frobenius manifolds [22]. 

5 Hodograph equations for hydrodynamic type sys- 
tems as equations for critical points and Euler- 
Poisson-Darboux equation 



Critical points for unfoldings of A and D singularities as the functions of deformation 
parameters represent very particular solutions of the BH equation and the system (64). 
What about the other solutions of these integrable equations? Do they describe certain 
unfoldings of critical points for other singularities? 

To address this question we begin with the BH equation = {u^)ti- A standard 
hodograph equation for it (see e.g. [36, 38]) 

h + 2t2U + f{u) = (68) 

is the equation defining critical point u{t) of the function W{x;t) = Fq{x) + t2X^ + tix 
where FQrc{x) = f{x) and f{u) is the function inverse to the initial data function Mo(ti) = 
u{ti, ^2 = 0). Let us consider a family of BH solutions which corresponds to the family of 
initial data with Fq{x) = x""*"^ + tnX"' + tn-ix"'~^ + ... H-tsX^ where t^, ...,tn are parameters 
and n is arbitrary. These solutions u{tit2, ts, tn) of BH equation clearly are associated 
with unfoldings ( versal at t,„ = 0) of An singularities considered in the previous section. 
Thus, the class of solutions of BH equation which corresponds to the family of initial 
data of this type with unbounded n describes unfoldings of critical points for all An 
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singularities. Alternatively, the parameters t^, ...,tn can be viewed as times for higher BH 
equations within the infinite BH hierarchy (n — )■ oo) [29]. 

An extension of this construction to the system (64) is almost straightforward. We 
already saw that the particular solution of this system describes motion of the critical 
points for the function (60). Thus, we will look for the functions F{x,y]t) = FQ{x,y) + 
t^{ax^ + y"^) + tix with Fq different from ax^ + Sxy"^. It is convenient to pass to the 
variables xi = x + y, X2 = x — y. In these variables F = Fo{xi,X2) + ^((^ + + 
2exiX2 + {e + l)x|) + ^i^^ + ^2) where £ = |(a — 1) and ^2 = 2^3. Equations defining the 
critical points /Sj^and are 

= 0, = 0. (69) 

Lemma. Solution of equations (69), with the function F of the form given above and 
any function Fq such that -F^g^/Sj = and Ff^^i^^Fj^^ii^ 7^ 0, is a solution of the system (67). 

Proof. Differentiating equation (69) w.r.t. ti and t2 and assuming that Fp^p^ = 0, one 
gets 

Pu, = -7T^, /32., = -777^ (70) 



Of? ^"-l O f? 



and 



_ ((e + l)/3, + e/3,) _ {ef3, + {e + 1)(3,) 

Eliminating -F/Sj/?^ and Fj^^/s^ from equations (71), one obtains the system (67). □ 



The function Fq due to the relations 



/3u,(ti,t2 = 0) = -— . , f3,,^{h,t2 = 0) 



^-ro/3i/3i 1*2=0 ^-r 0/32/32 1*2=0 

is defined implicitly by the initial data for and ^2- The functions F and aF where a is 
an arbitrary constant, obviously, give rise to the same equation. 
Equations (67) imply that 

1^1 + P2)t2 = \{i^ + + 2^/3i/32 + + . 

So, |(/3i + /32) = <^ij, i {{e + l)(3l + 2e(3^(32 + {e + l)/3^) = (^^^ where (p is a function of h 
and t2 and hence 

F{f3„f3,-dt) = d^. (72) 
There is a special subclass of functions F of particular interest. It is given by 

where 7 denotes a large positively oriented circle on the A plane. Written explicitly the 
function W is the series 
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W = -ti{xi+X2) + t2-[{e + l){xl + xl) + 2exiX2] + 

+ 1) [{e + 2)(x? + xl) + 3e{xixl + X2xl)] + (74) 
+t4^(e + 1) [{e + 2){e + 3)(xf + xt) + 4e(e + 2){xixl + xax?) + 6e{e + l)xlxl] + ... 

An important property of the function W is that it obeys the Euler-Poisson-Darboux 
equation E{e,e), i.e. 

(xi - X2)W^,^, = eiW,, - W,,). (75) 

This equation and representation of its solutions in the form (73) are known for more 
than a century (see [40]). In the papers [33, 34, 35] it was observed that the hodograph 
equations for the one-layer Benney hierarchy {e = \) and dispersionless dToda hierarchy 
{e = — |) are nothing but that the equations for critical points of the function 2eW. 
For arbitrary e equations 

= 0, Wp^ = (76) 

for the critical points are the hodograph equations for the system (67). Due to equation 
(75) the function W at 7^ automatically verifies the condition Wj3_^p^ = 0. For 
the function F of the form (73) the function Fq is a special one. The variables ^3,^4, ... 
can be viewed as the variables parametrizing a family of initial data for the system (67). 
Alternatively, one can treat them as the higher "times" for the commuting systems 

(3u, = 0ik{(3,,(32)(3u,, (32t, = 02ki(3„(32)(32u, A: = 3,4,... (77) 

with the characteristic velocities 9ik{(3^,(32) = ^{^)\x=i3, 02k{Pi,P2) = af;(S^)U=/3- 
The totality of these systems is an infinite hierarchy of systems associated with the system 
(67) and the function W plays the role of the generating function for this hierarchy. 
In the singular case e = the generating function W£=o is given by 




which gives rise to two independent BH hierarchies for and /32. 

The observation that the hodograph solutions of the two-component hydrodynamic 
systems of e-type describe critical points of function (73) is extendable to the multi- 
component case. It was shown in [33] that the critical points of the function 
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at £ = I are described by hodograph solutions of the dispersionless n-component coupled 
KdV equation. 

For arbitrary e 7^ the function (79) obeys the Euler-Poisson-Darboux system 

(xi = - W^xJ, i ^ k;i,k,= l,...,n. (80) 

From the equations Fp^ = 0,i = 1, ...,n for critical points (3i, of the function (79) 

one finds (3-^^ = = l,...,n;k = 1,2,... where Wk are defined by the expansion 

W{x\t) = J2k>i^kWk{x). Since = 1 one gets the following hierarchy of hydrody- 

namic type systems governing the motion of the critical points 

Pu, = Wkp.W)P,,^, ^ = l,...,n; k = 2,3,... (81) 

The first member of this hierarchy is given by the e-system 13^^^ = {£{J2m=i Pm) + 
f^i)f^iti^ ^ = 1,...,?T, considered in [39]. The system (81) implies that Wkti — ^it^ = 
0, k,l = 1, 2, .... So, Wk = hence 

WiP; dt) = dip. 

It was observed in [31] ( see also [39]) that the systems (81) have an interesting property: 
densities P of their conserved quantities satisfy the Euler-Poisson-Darboux equations dual 
to (80) ( i.e. with opposite sign of e). An infinite family of such densities is given by 

7 

A class of the hydrodynamic type systems for which hodograph equations coincide 
with the equations for critical points of certain functions is, in fact, larger. It consists of 
all semihamiltonian diagonal systems for which characteristic velocities 61 are components 
of the gradient of a function, i.e. 

/3;, = $^,/3;„ l = l,...,n. (82) 

where ^'/j, -^ Indeed, according to [11] the generalized hodograph equations for semi- 
hamiltonian diagonal system 13^ = 9i{f3)(3i^,l = 1, ...,n are given by the system 

x + ei{^)t + uji{/3) = 0, l = l,...,n (83) 
where the functions u!i{/3) obeys the equations 

— = a a ^ ^ ^ ^ (84) 

If 6i = $/3.then the l.h.s. of (84) is skew symmetric and hence Ui^^ = ujkp^- So there exists 
a function such that uji{[3) = ^uji3.,i = 1, As a result, the hodograph equations 

(83) take the form 

X + $^^t + = 0, l = l,...,n 
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that coincides with the equations Fg^ = for critical points of the function 

W = x{xi + ...Xn) +t^{x) + ^^{x). (85) 
From (84) it follows that this function obeys the system of equations 

W..., = ^^^^{W..-W.,), i^k- iX=l,...,n. (86) 

Xi, X}^ 

as well as the function ^i^. Densities P of the conserved quantities for the system (82) 
satisfy the equations [41] 

Pm. = - ^^''% (Pp.-Ph)^ ^^k- z,A.,= l,...,n. (87) 

Thus, the hydrodynamic type system (82) describes critical points of the function W (85) 
which obeys equations (86) and this equations (as well as the equations for the generating 
function of symmetries) and the equations for conserved densities P of (82) are dual 
to each other. 



6 Hierarchies of integrable systems and normal forms 
of singularities 

Now let us discuss the hierarchies of two- and three-component integrable systems con- 
sidered in the previous section and the corresponding functions W from the singularity 
theory viewpoint. The functions W = Ylkyi^f'^k provide us with an infinite families of 
symmetric homogeneous functions Wk{xi, ...,Xn) of degrees k. Is there any relation be- 
tween these functions and normal forms of germs in singularity theory? We will present 
here two observations which indicate that such a connection exists. 

In the two-component case the function W (73) in the variables x and y is given by 



W{x,y-t) = 




or 



W = tix+-t2{ax'^+y'^) + -{a+l)t3{ax^+3xy^)^ -{a+l)t4{ax^+bx^y'^+cy'^) + ... (89) 

2 6 384 

where a = 15a^ + 24a - 15, b = + 8a - 33, c = 6a^ - 48a - 102. 

Third term in (89) generates the first higher commuting flow for the system (64). At 
the same time, it provides us with the normal form of simple D4 singularity after the 
trivial rescaling of x. Fourth term gives rise to the next commuting flow for (64). Does it 
correspond to some standard normal form? It is easy to see that, rescaling x and y, one 
can convert it into the form x'^ + dx^y'^ + y^ where d is certain function of a . For generic 
a the parameter d 7^ ±2. Thus one has the normal form of the unimodular Xg singularity 
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(see [42, 4]). Higher Wk in the expansion W = X]fc>i^^'^'s which generate higher flows 

are of the form Wk = ^;=o c^z^;^"^'?/^' where ai{a) are certain polynomials in a. So, they 
describe unimodular singularities of the order k with the germs symmetric with respect 
to the reflection y — j- —y. 

Second example is given by the system (81) with n = 3 and e = ^, i.e. by the 
dispersionless three-component coupled KdV system. The function W is 



W = ■^ti(Xi + X2 + X3) + it2(3xj + + + 2XiX2 + 2x1X3 + 2x2X3) + 
2 

+—^3(5x1 + 5x2 + Sa^Q + 3x^X2 + 3x^X3 + 3X3X1 + 3X3X3 + 3X3X1 + 3X3X2 + 2x1X2X3) 
16 

In terms of variables x,y,z defined by 

xi = ax + y + z, X2 = ax + qy + q^z, X3 = ax + q'^y + qz, (90) 
where q=exp(^) and a = -3^, one has 

W = ^ahx + ^t2(5aV + 4yz) + ^h{x^ + y^ + + '^axyz) + ... (91) 

The transformation (90) is, in fact, the well-known relation between roots Xj of a cubic 
equation and its Lagrange resolvents ii modulo the identification ax = ^ii,y = \^2-,z = 
\h (see e.g. [13]). 

The term in (91) represents the normal form of the unimodular parabolic singu- 
larity Pg ( or T3 3 3 ) of corank three (see [42, 4]). At = 0, A; = 4,5,... the function 
(91) gives the unfolding of Ps singularity. The dependence of the critical points u,v,w for 
this function on deformation parameters ti and ^2 is described by the three-component 
dispersionless coupled KdV system (81) or equivalently by the system 

law, -w, -V 
aw, |m, W 
aw, V, |m 

The functions Wk with k = 4,5, ... represent higher order singularities of corank three 
and the function (91) give their unfoldings. 







) : 


(92) 
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